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We analyze the active-galaxy correlation reported in 2007 by the Pierre Auger Collabo- 
ration. The signal diminishes if the correlation-function approach (counting all "source- 
event" pairs and not only "nearest neighbours" ) is used, suggesting that the correlation 
may reveal individual sources and not their population. We analyze available data on 
physical conditions in these individual correlated sources and conclude that acceleration 
of protons to the observed energies is hardly possible in any of these galaxies, while heav- 
ier nuclei would be deflected by the Galactic magnetic field thus spoiling the correlation. 
Our results question the Auger interpretation of the reported anisotropy signal but do 
not contradict to its explanation with intermediate-mass nuclei accelerated in Cen A. 
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1. Introduction 

Astronomy of ultra-high-energy (UHE; energy E > 10^^ eV) cosmic rays (CRs) is a 
new branch of science whose birth we are currently witnessing. To identify UHECR 
sources is therefore an important and challenging task. Observational evidence in 
favour of particular sources is limited by low statistics, by systematic uncertainties 
in determination of the primary-particle type and energy, by deflection of charged 
particles in unknown cosmic magnetic fields and by other confusing effects in the 
analysis. Theoretically, several compelling mechanisms of UHECR acceleration have 
been developed and a number of potential astrophysical sources have been suggested 
(see e.g. Refs. for reviews and summary). Physical conditions in potential sources 
are often uncertain thus preventing one from firm identification of theoretically most 
favoured candidates. 



1 



30, 2010 2:28 WSPC/INSTRUCTION FILE sourcesIIv2 



2 S. Gureev and S. Troitsky 

Recently, considerable interest has been attracted by a claim of correlation of 
positions of nearby active galaxies with arriva l dir ections of UHECR particles de- 
tected by the Pierre Auger Observatory (PAO) The correlation was interpreted 
as an evidence that UHECR particles are protons either from nearby active galac- 
tic nuclei (AGN) or from other sources with a similar distribution in space. The 
AGN interpretation was considered likely based on the previous works favouring 
AGN as possible UHECR accelerators. Notably, Ref. SI cites Ref. ^ in this context; 
the latter reference however studies the most powerful AGN (BL Lac type objects 
and optically violent variable quasars), while the analysis of Ref. ^ is based on the 
nearby objects listed in the Veron-Cetti&; Veron catalog ^ which are mostly low- 
power Seyfert galaxies. In Ref. we analyzed general constraints from geometry of 
the source and from radiation losses which restrict potential astrophysical accelera- 
tors and compared them to the most recent astrophysical data. Our study suggested 
that only the most powerful AGN are able to accelerate UHE particles while Seyfert 
galaxies generically are not. In this work, we apply the general constraints discussed 
in Ref. to particular active galaxies correlated with Auger events. 

The rest of the paper is organised as follows. In Sec. [21 we br iefly review the Auger 
result emphasising particular details of the analysis of Refs. and demonstrate 
that the correlation signal is strong only if a single nearest-neighbour potential 
source is counted for each event and diminishes in the correlation-function approach 
when every pair "event-source" is counted both in the data and in the simulations. 
We then select these nearest-neighbour source candidates (mostly Seyfert galaxies) 
for the Auger events and study, in Sec. [3J their properties. We conclude that it is 
unlikely that these objects can accelerate protons up to the observed energies and 
discuss implications of this conclusion in Sec. SI 

2. Summary and discussion of the Auger correlations 
2.1. Anisotropy reported by PAcB 

Recently, the Pierre Auger Collaboration reported a significant deviation from 
isotropy of the arrival directions of the highest- energy cosmic rays observed by 
their surface detector 131. They studied an (unpublished) set of events with energies 
E > Ax 10^^ eV and cross-correlated their arrival directions with positions of active 
galaxies from the Veron catalog Three cuts were tuned to maximize the signal: 
(1) the minimal reconstructed energy of a cosmic-ray particle, £'inin; (2) the maximal 
redshift of a galaxy, Zmax! (3) the maximal angular separation ip between a galaxy 
and a cosmic-ray arrival direction at which they are counted as correlated. An excess 
of correlated event-galaxy pairs over the number expected for isotropic distribution 
was found; the signal was maximized for Smin = 5.6 x 10^^ eV, Zmax — 0.018 and 
iP = 3.1°. 

After the signal had been verified in an independent data set with the same cuts 

■^See Note Added for discussion of the most recent results. 
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(the probability of chance correlation in the independent set of 12 events was found 
to be P ~ 1.7 X 10"'^), details of the correlations were studied in the full data set HI. 
The full set was again subject to the same scans in -Emin, -Zmax and V'; the minimal 
formal probability of Pmin ~ 5 x 10~^ occured at Pmin = 5.7 x 10^^ eV, Zmax = 0.017 
and ip = 3.2°. The part of the data set which maximizes the correlation signal (27 
events) was published - and may be used for analysis and further studies while 
the full data set remains unavailable. To estimate P from Pmin, one may use the 
so-called penalty factor (see e.g. Ref. ^ which compensates for the artificial signal 
obtained by tuning the cuts or, in other words, for multiple tries (see also Ref. for 
discussion). This penalty factor should be determined by Monte-Carlo simulations^^ 
which are however quite lengthy for these low Pmin- One may estimate it roughly 
by counting the number of tries: the scan in < Zmax £ 0.024 in steps of 0.001 
gives 25 tries, the scan in 1° < V' < 8° in steps of 0.1° gives 71 tries and the scan in 
Emin < 4 x 10^^ eV with unspecified step gives effectively the number of tries equal 
to the number of events in the sample, that is 27. Altogether, the penalty factor is 
roughly 25 x 71 x 27 w 4.8 x 10"* (in practice, the penalty factor should be somewhat 
lower because the tries are not independent) and Pest ^ Pmin x 4.8 x 10* ?» 2 x lO^**. 

The consistency in probabilities (Pest '~ P up to an order-of-magnitude cor- 
rection which may be attributed to unpublishetiB tries with other catalogs) and in 
the values of the tuned cut parameters in the full data set and in the first half 
supports that a significant deviation from isotropy, and not an outstanding chance 
fluctuation, was indeed observed. The interpretation of this anisotropy is however 
not unique. Firstly, active galaxies may serve only as a tracer of the actual sources 
which could be distributed in the Universe in a similar way — indeed, most potential 
sources follow the distribution of galaxies which locally is not uniform due to super- 
clusters and voids. This interpretation cannot be distinguished from the simplest 
hypothesis that the AGN are cosmic-ray sources based on the present-level corre- 
lation signal; both are unified in the AGN hypothesis of the Auger collaboration: 
most of the cosmic rays reaching the Earth in that energy range are protons from 
nearby astrophysical sources, either AGN or other objects with a similar spatial 
distribution This interpretation requires a large (of order a hundred SI) number 
of sources. 

Secondly, it has been shown EE^HSUll that a similar signal may be easily pro- 
duced in case of just a few nearby sources but larger deflection angles (consistent 
with the chemical composition of the cos mic rays observed by the fluorescent detec- 
tor of the same Pierre Auger Observatory SSI and by the Yakutsk experiment In 
particular, the origin of a signifi can t part of cosmic-ray particles in the nearby FR I 
radio galaxy Cen A may explain the observed signal for the angular spread ~ 10° 
typical for intermediate- mass nuclei arriving from a source with these Galactic coor- 
dinates^^. Moreover, this kind of interpretation is actually supported by statistical 

''Or published elsewhere, e.g. in Ref. H^l. 
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analyses of the global distribution of the arrival directions of the cosmic rays in the 
maximally correlated Auger sample "'^1. Qualitatively, under the AGN hypothesis, a 
similar number of events is expected from the directions of the Virgo and Centaurus 
superclusters (the latter is significantly farther but the former is observed by Auger 
with much lower exposure). In practice, however, iio events came from Virgo while 
the Centaurus region dominates the correlation Alternatiye irrterpre of 
the Auger signal have been considered also in Refs. | 17 | 18 | 19 | 20 | 21 | 22 | 23 ] 

Thirdly, independent experiments disagree about the presence of a similar corre- 
lation in their data. Using a dataset of a similar size, the High Resolution Fly's Eye 
experiment (HiRes) did not find a correlation with nearby Veron AGN either by 
applying the PAO prescription or performing similar scans in i^minj -Zmax and ip'^. 
The data of the Yakutsk Extensive Air-Shower Array support the correlation 
the strongest signal occurs at iS'min, ^max and ijj close to the values prescribed by 
PAO though the angular resolution of the experiment and possibly the energy scale 
are different. We note that UHECR correlations with nearby Seyfert galaxies from 
catalog have been discovered by Uryson in 1999 "^^ using the AGASA and 
Yakutsk data for somewhat lower energies E > A x 10^^ eV. However, a similar 
sample of Seyfert galaxies from the Veron catalog does not correlate with the same 
cosmic-ray sample 

In this paper, we analyse the observed correlation from a completely different 
point of view. 



2.2. Nearest neighbour or correlation function? 

In any correlation analysis similar to that of the number of observed pairs "event- 
source" is compared to that expected from chance coincidence. The definition of the 
"number of pairs" is however not unique: when an event is found within the angular 
distance "0 from k sources and fc > 1, the "number of pairs" may be counted either 
as 1 (counting only nearest-neighbour source for each event) or as k (counting 
every source— event pair). Clearly, the same rules should be applied both to the 
data set and to the Monte-Carlo sets used to estimate the expected number. While 
the nearest-neighbour approach may be psychologically favoured because a single 
cosmic-ray particle cannot originate in several sources, it is the second approach 
which is widely applied in statistical studies by terms of the correlation function 
between two sets of points (events and sources 

For an isotropic distribution of potential sources on the sky (when clustering of 
sources at the angular scale ^ ip is consistent with chance coincidence), both ap- 
proaches give qualitatively similar results and both may serve to test the association 
between the ensemble of sources and the ensemble of cosmic rays. When sources 
are clustered at the angular scale ^ ip (which is the case for the AGN sample used 

•^This correla tion-fun ction approach has been used in many previous UHECR correlation studies, 
see e.g. Refs. EESEQ] 
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by PAO), it is the correlation- function approach which should be used to study the 
correlation between ensembles, while the nearest-neighbour method may still be 
used to determine particular sources. If the correlation exists with the population 
of (clustered) sources, one cannot say which particular objects are responsible for 
individual events; the correlation is a collective effect of many sources separated 
from the cosmic-ray arrival directions by angles ^ ip. A similar effect is expected if 
the sample of the sources under study is just a tracer of other sources with a similar 
distribution. On the other hand, if the correlation-function approach reveals no sig- 
nal but the nearest-neighbour approach does reveal some, this is a strong argument 
that some of these nearest-neighbour sources may be related to particular events 
(we will see below that it is precisely what happens with the Auger correlation). 
This statement, to be discussed in detail elsewhere, may be easily understood by 
comparison of two cases: a large number of weak sources versus a small number of 
strong sources in the sample. 

The key poin t of the calculation of the probability of chance correlation used in 
PAO papers 1^1^ is the calculation of the fraction of the (exposure- weighted) sky area 
covered by circles with radius ^ cen tered on the potential sources in the catalog. 
This fraction, called p in Refs. allows to estimate the probability of chance 
correlation in the nearest-neighbour approach simply by the binomial distribution, 
cf. Ref. As it is illustrated in Fig. [1] p may be calculated in two different ways: 
in the "nearest-neighbour" approach, the total covered area of the sky is counted 
independently of overlaps, while in the correlation-function approach, the overlap 
area of two circles is counted twice. 

Let us study, by means of the two methods, the sets of 27 cosmic-ray events 
and of 422 AGN which maximize the Auger correlation. In the nearest-neighbour 
approach and for tp = 3.2°, we find p w 0.21, the number of correlated pairs is 20 at 
5.6 expected by chance, so Pmin ~ 5 x 10^^ in agreement with the Auger claim. In 
the correlation-function approach, p w 0.336, the number of correlated pairs is 25 at 
8.4 expected by chance, so Pmin ~ 1.7x 10"'^, several orders of magnitude highei0. In 
the case under discussion, these two differ qualitatively: when Pmin is multiplied by 
the penalty factor ^ W^, the nearest-neighbour mehtod gives, as discussed above, 
a significant signal, while the correlation-function result is consistent with absence 
of any correlation. 

This particularly strange pattern of correlation signals suggests that, if the signal 
is not a chance fluctuation, then the "nearest neighbours" of correlated events should 
be the actual sources of cosmic rays. In the rest of the paper, we will therefore study 
the properties of these "nearest neighbours" . 



"^In the latter case, the binomial formula used in Refs.EElig not applicable; Pmin is estimated by 
the Monte-Carlo simulations. 
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"Nearest neighbour" "Correlation function" 




S<3So S=3So 
(1 pair) (3 pairs) 



Fig. 1. Illustration of two approaches to estimate the correlation signal for clustered sources. 
Suppose that the circles of radius i/i (area So) around fc > 1 potential sources (grey circles on the 
plot; fc = 3 in this example) overlap and the cosmic-ray arrival direction (star) falls in the over- 
lapping region. In the nearest-neighbour approach (left panel), each cosmic-ray event is associated 
with only one source; the area S covered by circles on the sky is counted once independently of 
overlaps. In the correlation-function approach (right panel), each pair "source— event" is counted 
(fc = 3 pairs in the example) and the overlapping area is counted fc times in the calculation of S. 

3. Physical conditions in the correlated galaixies 

The list of the sources responsible for the Auger correlation is given in Table [TJ In 
this section, we analyze available data on physical conditions in these galaxies and 
compare them with requirements for UHECR acceleration discussed in Ref. 

3.1. The central black holes 

We start with the analysis of conditions close to the central black holes of the 
correlated galaxies. Following the results of Ref. we estimate the maximal possible 
energy of UHECR particles with atomic mass A and charge Z from the central black 
hole of mass Mbh: 

For a few objects, the mass of the black hole has been estimated previously in 
the literature; for others, we estimate A/bh from available data. 

Methods of determination of the mass A/bh of the AGN central black hole have 
been reviewed in Ref. (for a recent update, see e.g.^^^, where basic corresponding 
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references can be found. For 3 of 17 galaxies correlated with Auger events, Mbh 
has been estimated in the literature using observations of circumnuclear dynamics 
by means of: 

(GD): resolved gas dynamics; 

(MM): observations of accretion-disk megamasers; 

(Ha): width of the Ha line. 

For other galaxies in the sample we used, depending on the data availability, the 
following methods. 

(SVD): stellar velocity dispersion. This quantity, a, was found to be corre- 
lated with Mbh; we use the relation ^^ll 

For some of the galaxies, a is given in the LEDA database'^; for a few others, we 
found it in the literature. 

(KSB): buldge Kg magnitude. It is also correlated with A/bh; we use the 
relation 

log (y^^) = (-2-5 ± 0.6) + (-0.45 ± 0.03) 

where the absolute bulge magnitude K'^ is related to the observed bulge magnitude 
Kf through 

X° = + Ak + DM. 

Here Ak is the correction for the Galactic extinction in the K band I^HI and DM 
is the distance modulus (both Ak and DM are given in NED). For some of the 
galaxies, was determined by the disk-bulge decomposition in previous studies; 
when known, it is often quoted in NED. 

(KSM): bulge Kg from morphology. When is not directly available, one 
may estimate it from the total Kg magnitud e a nd the Hubble morphological type 
T of a galaxy by using the following relation ^35]^ 

= Ks+ 0.297 (T + 5) - 0.040 (T + bf + 0.0035 (T + 5f . 

The precision of this relation is of order 0.5™ and depends on T '^Sj. While this 
method is the least precise, it may be used to estimate Mbh for any galaxy with 
known distance because either or a lower limit on it may be determined from 
the 2MASS full-sky survey data. For 9 of 17 correlated galaxies, KSM is the only 
available method to estimate Mbh- 

The data we use are collected in Table [2] while the summary of A/bh estimates 
is given in Table [3l For estimates of the maximal energy of cosmic-ray particles we 
always use the most precise estimate of the black-hole mass given in the bold face 
in Table |3l When different methods are applicable, less precise estimates are in a 
good agreement with better ones (a comparison plot is given in Fig. [2|). 

Figure [3] compares the maximal attainable energy -Emax tor protons and nuclei 
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KSM SVD SVD Hfi- 



KSM GD SVD GD 



Fig. 2. Different estimates of the mass of the central black hole in two galaxies studied in Sec. [3] 
Letters denote the method of estimation (see Sec. I3.l(l . 




Log(Eobs/eV) 



Fig. 3. The maximal energy -Emax attainable by protons and carbon nuclei in the vicinity of 
the central black hole of individual active galaxies correlated with the Auger events, versus the 
energy -Eobs of the observed cosmic-ray particle associated with the same galaxy (if a source 
corresponds to two events, -Eobs is the highest one). The shaded area corresponds to the allowed 
region i?obs ^ ^max for protons; the area above the thick black line corresponds to the allowed 
region i?obs £'max for carbon nuclei. The error bars correspond to precision of determination 
of the black-hole mass; Znajeck conservative upper limit for the magnetic field is assumed (see 
Ref. Efor details). For PC 2207-1-0122, a dwarf galaxy associated with two cosmic-ray events but 
not detected in the 2MASS survey, only an upper limit on i?max is available. 

near the central black hole of the galaxies with actual energies of the correlated 
cosmic-ray events. The energy i^max has been calculated using Eq. ([T]). We remind 
that this equation gives an upper limit on E'max which is lower in particular models 
and for realistic values of the magnetic fields. Central black holes of any of these 
galaxies cannot accelerate protons to the observed energies. While they could be 
able to accelerate heavier nuclei, deflections of these nuclei in the Galactic magnetic 
field would spoil the directional correlation at ■0 3°. 

3.2. Jets and other extended structures 

In eight of the correlated galaxies, extended nuclear outflows have been detected. 
Except for Cen A, which is discussed separately in the following section, these 
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outflows are most probably non-relativistic, and when the origin of the emission 
may be determined, it is thermal. Observations of outflows in correlated Seyfert 
galaxies are summarized in Table |4l Magnetic fields in these outflows have not been 
studied; to get order-of-magnitude estimates, we assume that the energy of the 
magnetic field is the same as the total energy of the radiation field. We assume a 
spherical source of the radius R/2 filled with magnetic field i?, so that the energy 
of the magnetic field is Ei ^ ttR^B^ /&. The energy of the radiation field may be 
estimated as E2 ~ RdE2/dt, where the energy loss dE2/dt = AttD'^ jl^dv for a 
source at the distance D with the observed flux at the frequency v (we use the 
definition of Ii, as the energy per unit area per unit time per unit frequency interval 
per unit solid angle, that is is measured e.g. in Jansky). Within our precision, 
Jiudv ~ I^v. Equating Ei = E2 and expressing D/R through the angular size of 
the source 0, we arrive at 

1/2 

B - 6.5 X 10~« G 



mas / V'^Jy 

Estimated magnetic flelds (given in TableE]) do not exceed those expected for normal 
galaxies of the same size; the outflows certainly cannot serve as acceleration sites 
of UHECRs because they do not satisfy the Hillas criterion. 

3.3. The case of Centaurus A 

One of the correlated galaxies, NGC 5128, hosts a radio galaxy known as Cen A. Due 
to its proximity, it represents a textbook example of FR I radio galaxy, st udied in 
great deta ils for decades. It has been suggested long ago as a source of some l^^^^ or 
even most^^H^of the obs erved cosmic ra ys of extreme energies. Both the originally 
and some of the alternativel '*^'*' ! '*^^ ! '*''^ ! '*^^ ! '*^^^ '^'^ interpretations of the Auger correlation 
signal imply that Cen A is a UHECR source. We therefore devote a separate section 
to the study of physical conditions in Ce n A. 

Centaurus A (for reviews, see Refs. I ^'^ l ^^ l g^j^^j continuously updated web- 
page is a low-power radio galaxy (its luminosity is intermediate between 
Seyfert and FR I '^^^ which, most probably, experienced a major merging event 
recently. High-resolution studies determine the compact nucleus (optically classi- 
fled as Seyfert 2) whose multifrequency spectr al energy dis tribution is very similar 
to 



id as Seyfert 2) whose multifrequency spectr al energy dis tribution is very similar 
those of the BL Lac type objects (see Refs. I 56 | 57 | 58 [59j discussioiR) ; the jet 



which experience several bendings observed at scales from light months'^ to dozens 
of kiloparsecs from the nuclei; inner radio lobes associated with the jet and giant 
outer radio lobes which span about 9° on the sky. The jet is mildly relativistic, 
r < 2.5 at parsec scales (velocity ~ 0.45c was inferred and proper motions 
^ 0.12c were observed 1^^. Polarization studies reveal that the magnetic held in 



°Cen A is the onlvpbject in the AGN sample studied by PAO which is classified as a BL Lac in 
the Veron catalog!^. 
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Log(size/kpc) 



Fig. 4. The size-field plot for different parts of Cen A. Grey colours: the nuclear region (light- 
grey er ror- bar diagonal, radio-to- X-ray modelling of Ref. ESI dark-grey error-bar diagonal, that 
of Ref. ES grey box, synchrotron self-absorption measurement of Ref. ESI). Gre en colours: the 
jet (small green triangles, minimal-pressure estimates for resolv ed components ESI light-green 
empty tria ngle s, equipartition estimates for resolved components ED green box, model esti mat e 
for a knot 1221 green diamond, model estimate for the possible hard-spectrum shear layer ESI). 

Blue col ours: the lobes (large light-blue triangle, equipartition estimate for the giant outer 

lobes EH (see also Ref . ESI) ■ other points arc estimates for inner lobes: small light- blue triangle, 
minimal-energyESl small empty light-blue triangle, equipartition dark-blue triangles, minimal- 
pressure ESI dark-blue empty triangles, equipartition ED) . The allowed region for acceleration of 
7 X lO'^^ eV particles is located between red lines (thick, protons; thin, carbon nuclei; dashed, iron 
nuclei; lower lines corresponds to the Hillas limit, upper ones corresponds to the radiation-loss 
limit for diffusive acceleration). Note that the diffusive-acceleration loss limit is determined by 
A/Z and is therefore indistinguishable for iron {A = 56, Z = 26) and carbon [A = 12, Z = 6). See 
Ref. 1^ for details of the calculation of allowed areas and for description of the methods to estimate 
the magnetic fields. 



the jet is parallel to the jet axis at the scales of at least 3 kpc (see e.g. Ref. 
and references therein). Several dosen knots are resolved in the jet; their radio and 
X-ray positions are systematically offset from each other (see e.g. Ref. These 
knot s ar e possible sites of particle acceleration; however recent X-ray observations 
(e.g. suggest that another, distributed along the jet, acceleration mechanism 
should be at work. This mechanism may be associate d wi th a hard-spectrum shear 
layer possibly observed along the kiloparsec-scale jet Multifrequency observa- 
tions (e.g. Ref. establish the synchrotron origin of the jet emission with high 
confidence. 

Numerous estimates of the magnetic field in different parts of Cen A are sum- 
marized in Fig. m Acceleration of protons to the energy of 7 x IQ-'^^ eV (the energy 
of the events associated with the Cen A nucleus by PAO) is hardly possible in this 
source, except maybe for the giant outer lobes which are however displaced from 
the nucleus by an angular distance exceeding ip = 3.2°. On the other hand, acceler- 
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ation of intermediate-mass and heavy nuclei in the jet and lobes of Cen A remains 
a plausible interpretation of the correlation signal. 

4. Conclusions 

The correlations observed by PAO may be explained either by the AGN hypothesis 
or by chance coincidence. The AGN hypothesis implies that the primary cosmic-ray 
particles are protons and has three options: 

(i) the correlated AGN are (some of) the sources; 

(ii) the AGN as a whole are sources; particular correlated AGN may not be the 
actual sources if they are located in regions of enhanced AGN density; 

(iii) the sources are not AGN but follow a similar distribution in space, that is, 
like AGN, they follow the local large-scale structure of the Universe. 

The options (iii) and (ii) are dis fav oured by the absence of the events from the 
direction of the Virgo supercluster El the options (i) and (ii) are disfavoured by 
the luminosity argument ^1^. 

In this paper we have noted that (ii) and (iii) are also disfavoured by the dis- 
appearence of correlations in terms of the correlation-function approach, contrary 
to the nearest-neighbour probability estimate. In order to test the option (i), we 
have studied physical conditions in particular correlated galaxies. Most of them are 
low-power Seyfert galaxies; we have demonstrated that none of the correlated 
AGN can accelerate protons to the observed energies, thus excluding the 
option (i) and giving further support to the explanation of correlations by means 
different than the AGN hypothesis. 

Our results suggest that the correlated galaxies, notably Cen A, can accelerate 
nuclei to the observed energies. The deflections would however spoil the small-angle 
correlations, but the origin of a significant part of the detected particles in Cen A jet 
and lobes remains one of the most probable explanations for the observed anisotropy. 

Note added 

Recently, a few month s af ter this paper was completed and posted in arXiv, PAO 



new data. In agreement with our results, these data speak against the AGN in- 
terpretation of the observed anisotropy: among 31 events with E > i?min recorded 
after publication of Ref. ^ and before March 2009, 8 correlate with AGN while 6.5 
chance coincidences are expected for isotropic distribution of arrival directions 
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Table 1. Active galaxies correlated with Auger events. Col. (1) gives the source name, 
Col. (2) gives the energy(ies) of the cosmic-ray particles for which this source is the "nearest 
neigh bou r" in the correlated sample; Col. (3) gives the distance to the object (taken from 
NED EH); in Cols. (4) to (6) objects detected in radio (R), X rays (X) and gamma rays 
with E > 1 MeV (7) are marked with +; Col. (7) gives the AGN type (mostly from NED: 
Seyfert galaxies (Sy), galaxies with H II regions (HIl), emission- line dwarf galaxy (EmG) 
and FR 1 radio galaxy (FRI)); Col. (8) gives the estimated mass of the central black hole 
(see Sec. I3.ll and Table |3] for details and error bars). Objects with relativistic jets (J) and 
non-relativistic outflows (O) are marked in Col. (9), and the object with knots (K), hot 
spots (H) and lobes (L) is marked in Col. (10). 



Object 


E 


d 




detection 




AGN 




jets 


knots. 




(EeV) 


(Mpc) 


R 


X 


7 


type 




lobes 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


ESO 383-G18 


84 


51 


+ 


+ 




Sy2 


6.60 







4U 1344-60 


66 


51 




+ 




Syl.5 


6.39 






ESO 139-G12 


83, 59 


68 








Sy2 


7.23 






IC 4518A 


63 


66 


+ 


+ 




Sy2 


5.77 






NGC 424 


84 


46 


+ 


+ 




Sy2 


7.34 







NGC 4945 


58 


5 


+ 


+ 




Sy2 


6.15 







IC 5169 


57 


42 


+ 






Sy2 


7.52 







CGCG 374-029 


59 


58 


-1- 






Syl 


6.01 






NGC 1346 


85 


54 


-1- 






Syl 


7.38 






NGC 7591 


83 


69 


-1- 






Sy 


7.58 






NGC 1358 


69 


54 


-1- 


+ 




Sy2 


7.71 







Cen A 


69, 70 


5 


-1- 


+ 


+ 


FRI 


7.65 


J 


K,H,L 


PC 2207-1-0122 


58, 71 


55 








EmG 


<5.76 






NGC 2989 


64 


56 


-1- 






HII 


6.79 






NGC 1204 


78 


57 


-1- 






HII 


7.75 






NGC 7130 


90 


65 


-1- 


+ 




Sy2 


7.35 







NGC 5506 


83 


29 


-1- 


+ 




Sy2 


6.70 
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Table 2. Data used to estimate the mass of the central black hole in active galaxies 
correlated with Auger events. Col. (1); the source name. Col. (2): Kg magnitude (A: 
95% confidence-level lower limit from t he 2MASS point source rejection table; B: 
bul ge d ecomposed magnitude from Ref. EZ| C: buldge decomposed magnitude from 
ReflSl, D: DENIS value from LEDA, no letter: 2MASS value from NED; for B and C, 
2MA SS e rror bars assumed). Col. (3): correction for Galactic extinction in K band 
(Ref. EH, NED). Col. (4): distance modulus (NED; accuracy is ±0.15"" except for 
PC 2207+0122 (±0.21™). Col. (5): Hubble morphological type from LEDA. Col. (6): 
central velocity dispersion (star denotes data from Ref. for which 10% accuracy 
was assumed in calculations, if no reference given, data from LEDA). 



Object 


Ks, mag 


Ak, mag 


DM, mag 


T 


o", km/s 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


ESQ 383-G18 


11.99 + 0.05 B 


0.022 


33.52 





92.3 * 


4U 1344-60 


11.140 ± 0.043 


1.082 


33.49 


5 




ESQ 139-G12 


9.709 + 0.039 


0.027 


34.17 


4 




IC 4518A 


11.09 + 0.22 D 


0.058 


34.10 


6 




NGC 424 


11.23 + 0.02 B 


0.006 


33.31 





142.6 * 


NGC 4945 


4.483 + 0.017 


0.065 


28.43 


4 


127.9 + 19.1 


IC 5169 


9.776 + 0.029 


0.006 


33.10 







CGCG 374-029 


11.260 + 0.050 


0.043 


33.80 


5 




NGC 1346 


9.790 + 0.029 


0.019 


33.66 


3 


116.1 + 4.7E3] 


NGC 7591 


9.602 + 0.032 


0.038 


34.19 


3 




NGC 1358 


10.18 + 0.003 C 


0.023 


33.65 





176.7 + 10.1 


Gen A 


3.942 + 0.016 


0.042 


28.39 


-2 


119.8 + 7.1 


PC 2207+0122 


> 17.057 A 


0.013 


33.69 


7 




NGC 2989 


10.229 + 0.049 


0.023 


33.72 


4 




NGC 1204 


9.947 + 0.030 


0.027 


33.77 







NGC 7130 


10.18 + 0.01 B 


0.011 


34.08 


1 


143.2 * 


NGC 5506 


10.74 + 0.015 B 


0.022 


32.29 


1 


97.9 * 
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Table 3. The estimates of the mass of the central black 
hole of active galaxies correlated with Auger events. Col. (1) 
gives the source name, Col. (2) gives the mass estimate. 
Col. (3) gives the method (see Sec. I3.l| l. Col. 4 gives the 
reference to either the mass value (v) or the data used to 
calculate it (d). If no reference is given, the LEDA value of rr 
was used in the SVD method, and the data of Table[2]were 
used in the KSB, KSM method. The most precise estimates 
for each object, which we use in this work, are given in the 
bold face. 



Object 


log ( Wlf-) 


Method 


Reference 










ESQ 383-G18 


- 1.40 ± 0.19 


SVD 






-0.80 ±0.88 


KSB 




4U 1344-60 


-1.61 ±0.92 


KSM 




ESQ 139-G12 


-0.77 ±0.94 


KSM 




IC 4518A 


-2.23 ±0.86 


KSM 


EZld 


NGC 424 


0.66 ± 0.18 


SVD 




-0.56 ±0.90 


KSB 




NGC 4945 


-1.85 ±0.18 


MM 






-0.98 ±0.93 


KSM 






-0.84 ±0.26 


SVD 






< -2 


GD 


[39]^ 


IC 5169 


-0.48 ±0.92 


KSM 




CGCG 374-029 


1.99 ±0.91 


KSM 


Ely 


NGC 1346 


-0.62 ±0.07 


Hq 




-0.75 ±0.90 


KSM 






-1.27 ±0.20 


SVD 






-1.01 ±0.09 


SVD 


ESld 


NGC 7591 


0.42 ±0.92 


KSM 




NGC 1358 


-0.29 ±0.10 


SVD 






0.14 ±0.95 


KSM 






0.07 ± 0.93 


KSB 


ESId 


Cen A 


-0.35 ±0.05 


GD 


HqIv 




0.22 ± 0.97 


KSM 






-0.95 ±0.12 


SVD 








GD 


Elv 


PC 2207-1-0122 


< -2.22 


KSM 




NGC 2989 


-1.21 ± 0.92 


KSM 




NGC 1204 


-0.25 ±0.93 


KSM 




NGC 7130 


-0.65 ±0.18 


SVD 


EZld 




0.26 ±0.93 


KSB 




NGC 5506 


-1.30 ±0.19 


SVD 


EZId 




-0.79 ±0.89 


KSB 






-0.25 ±0.20 


SVD 
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Table 4. Observations of outflows detected in correlated galaxies. Col. (1): the source 
name. Col. (2): linear size. Col. (3): frequency of observations. Col. (4): observed flux. 
Col. (5): estiroate of the magnetic field (for kpc wind in NGC 5506, from thermal source 
model of Ref. 33 for other cases, from Eq. lO). Col. (6): reference to the data. 



Object 


R, kpc 


u, GHz 


F, mJy 


i3cst, G 


Ref. 


(1) 


(2) 


(3) 


(4) 


(5) 




(6) 


ESO 383-G18 


0.107 


8.3 


1.7 


2 X IQ- 


-7 


\M 


NGC 424 


0.33 


8.3 


13.1 


10-7 




\M 




0.5 


1.5 


23.9 


5 X 10- 


-8 


\M 






8.4 


13.1 


9 X 10" 


-8 


Ha 


NGC 4945 


31.4 


4.8 


51 


9 X 10- 


10 


\M 




0.5 


1.5 X lO^** 


2.6 X 10"^ 


10-7 




\M 


IC 5169 


0.35 


1.5 


17.6 


5 X 10- 


-8 


Ha 






8.4 


3.7 


6 X 10" 


-8 


Ha 


NGC 1358 


0.36 


1.5 


3.4 


3 X 10" 


-8 


Ha 






8.4 


0.9 


3 X 10" 


-8 


Ha 


Cen A 






see Sec. 13.31 






H3 


NGC 7130 


0.188 


8.4 


18.1 


4 X 10" 


-7 


NGC 5506 


0.0015 


8.4 


24.4 


3 X 10- 


-5 


Ha 




0.36 


8.4 


67.6 


10-7 




H3 




1.2 ^ 3.0 


1.5 X lO^* 


(0.6-^ 2.4) X 10-7 


(0.8 -f- 6) X 


10-'^ 


H21 




0.302 


8.3 


95.0 


2 X 10" 


-7 





